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Abstract 

In the framework of the supersymmetric standard model, the lighter stau often 
becomes long-lived. Such longevity of the stau is realized in three well-motivated 
scenarios: (A) the stau is the next-to-lightest supersymmetric particle (NLSP) and 
the gravitino is the lightest supersymmetric particle (LSP), (B) the stau is the LSP 
and R-parity is slightly violated, and (C) the stau is NLSP, the neutralino is the 
LSP, and the their masses are degenerate. We study the event topology and the 
decay of the stopping stau at the hadron calorimeter at the LHC, and show that it 
is possible to identify the reason why the stau becomes long-lived. 



1 Introduction 



A supersymmetric (SUSY) standard model (SSM) is the most promising candidate of 
physics beyond the standard model. In some classes of the SSM, the lighter stau (fx) is 
long-lived. Such long-lived staus have great impact on both cosmology and LHC physics. 
Even if such a model is realized, the stau cannot be completely stable but must have a 
finite lifetime [TJ. If the decay length of the stau is longer than the detector size, the 
long-lived massive charged particle will be observed at the LHC [2j [3j HJ EE [6j. In this 
paper, we define the long-lived stau as the stau with the decay length longer than the 
detector size. By using the long-lived stau, it is possible to measure the property of the 
stau, such as the mass, lifetime and so on [5j [6j [7J [8j HJ [ID]. It is also possible to measure 
various quantities of the other particles using the long-lived stau [UJ [121 [13 ttH EE [TBI 

After discovering the stau signals, it is important how we probe the fundamental 
physics, using such observation. To do so, the most important thing is to identify the 
mechanism which makes the stau meta-stable. In general, a particle can be (meta-) stable 
because of symmetry, weak interaction and/or kinematical reason. In the framework of 
the minimal supersymmetric standard model (MSSM) with the gravitino, there are three 
well-motivated reasons why the stau can be long-lived: (A) the stau is the next-to-lightest 
supersymmetric particle (NLSP) and the gravitino is the lightest supersymmetric particle 
(LSP), (B) the stau is the LSP and R-parity is slightly violated, and (C) the stau is NLSP, 
the neutralino is the LSP, and the their masses are degenerate. Without knowing what 
makes the stau meta-stable, it may be impossible to probe the fundamental physics. 

In this paper, we discuss the discrimination of the three reasons by studying missing 
energy + stau SUSY events and the energy deposit from the stau decay which stops at 
the hadron calorimeter. We firstly describe the mechanisms which make the lightest stau 
stable in section [2j In section [31 we show that we can distinguish the case where the 
stau-longevity attributes to the kinematical reason from the other two cases by using the 
event topology. In section HJ we see that the other two scenarios, the gravitino LSP with 
the stau next lightest supersymmetric particle (NLSP) and the stau LSP with slightly 
broken R-parity can be also discriminated by studying the decay of the stau which stops 
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at the hadronic calorimeter. The last section is devoted to conclusions and discussions. 



2 Mechanism of Longevity of Stau 

The lighter stau can be long-lived by several reasons. Particularly, there are three well- 
motivated scenarios. 

(A) Weakly interacting LSP: 

When there is a weakly interacting LSP and the stau is the NLSP, the stau can be 
long-lived. For example, the gravitino and the axino are candidates of such an LSP. 
In this paper we discuss the case of gravitino-LSP. The lifetime of f± depends on the 
gravitino mass 7713/2, 



T, 



iM >t m k = 6 x 10 - sec (_!£_)- (^-Y , (1) 
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where Mp is the reduced Planck mass. If the gravitino mass is larger than 0(1) keV, 
T\ will be recognized as a long-lived charged particle at the LHC. 

(B) R-parity violation (RPV): 

If R-parity is violated, the LSP does not need to be neutral as long as it decays 
into the standard model sector with a finite lifetime. Hence, fi-LSP is possible. 
The lifetime of f\ depends on the size of RPV couplings, and f\ is long-lived if the 
couplings are small enough, fx becomes the LSP over the large parameter space of 
the gravity mediation scenario, and there exist allowed regions where fi-LSP with 
small RPV do not contradict to the experiments such as the observation of proton 
decay [30] and the cosmo logical bounds [3T] . 

(C) Kinematical suppression: 

Even if the R-parity is conserved and the lightest neutralino is the LSP as often 
assumed in the MSSM, the stau can be long-lived; the decay rate of the stau is 
highly suppressed if the mass difference between f\ and Xi is much less than the 
tauon mass. It has been discussed that such a scenario is favored in the context of 
solving the 7 Li problem |32j . 
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Depending on the lifetime of fi, the strategy to probe the origin of the stau longevity 
is different. If the typical decay length of f\ is longer than the detector size, f\ behaves 
as if it were stable in the LHC experiments, which is the case we study in this paper. In 
the case that the most of fi decay in the detector in their flight, a study of f\ kink track 
becomes more important. Such a study is found in the recently appeared paper [29J. 

3 LHC Event Topology in Long-lived Stau Scenarios 

Now we start our discussion to probe the origin of the stau longevity at the LHC. We 
concentrate on the case where the typical decay length of f\ is longer than the size of 
the detector system at the LHC, and in such a case the stau track can be detected by 
the detector. Since the stau is a heavy particle, its velocity is relatively small. The 
ATLAS detector and the CMS detector can detect such an exotic charged particle (2j [39] , 
relying on the fact that the most of the staus have small velocity. An existence of the 
slow-moving charged track with high-pr is the most characteristic signature of the SUSY 
event, and standard model backgrounds will be considerably suppressed by requiring the 
slow- moving track to the signal event. A part of staus still have a large velocity which 
may be almost the speed of light. Such a high-speed stau would not be distinguished 
from a muon. We refer the stau track with small velocity to the stau-like track, while the 
stau-track with large velocity and the ordinary muon track to the muon-like track. The 
SUSY events can be categorized by the number of fi-like tracks: (1) only one fi-like track 
is found, (2) two fi-like tracks are found. This categorization is useful in the following 
discussions. We do not consider the third category, no stau-like tracks in the event. 

First, consider a scenario that heavy superparticles, such as squarks, always decay 
into f\ via cascade decay. This is the case realized in, e.g., the scenario (A). Every 
SUSY event contains two staus from heavy superparticle decay, however, we may find 
only one stau-like track because we could not distinct the stau-track with large velocity 
from muons. Although f\ mimics a muon when its velocity is so large, such a fake muon 
generally possesses high-p^, ~ 0(100) GeV. Therefore, a high-pr muon will associate to 
the event with one fi-like track in this scenario. A similar discussion would hold for the 
scenario (B). If the RPV interactions are weak, heavy superparticles do not decay into 
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the standard model sector, and SUSY cascade decay chains will end up at the LSP, f\. 

In the scenario (C), however, the above discussion does not hold. In this case, heavy 
superparticles decay into not only fx but Xi with some fraction. When two SUSY cascade 
decay chains result in one f\ and one Xii we will detect only one fi-like track and, at the 
same time, observe a large missing transverse energy. In addition, it is expected that a 
small number of high-p^ muons are associated to the event with one fx-like track in this 
scenario. These signatures are different from the above scenarios, and hence they could 
be used to probe the origin of the longevity of the stau. 

In order to confirm the discussions, we consider three representative scenarios here: 
(A) fi-NLSP in mGMSB [33], (B) fi-LSP in mSUGRA with RPV, (C) mSUGRA with 
degenerated f\ and Xi masses. A mass spectrum of the scenario (A) is characterized 
by free parameters: {A, M mess , N 5 , tan/3, sign(/i)}. The scenario (B) and (C) have 
parameters: {m , M 1 / 2 , A , tan/3, sign(yu)}Jj We randomly generate parameter points for 
each scenarios over the following parameter spaces uniformly: 

(A) : {10 TeV < A < 500 TeV, A < M mess < 5000 TeV, 

1 < N 5 < 5, 1 < tan/3 < 60, sign(/i) = ±} (2) 
(B)(C) : {-3000 GeV < m < 3000 GeV, GeV < M 1/2 < 3000 GeV, 

-3000 GeV <A < 3000 GeV, 1 < tan/3 < 60, sign(/i) = ±} (3) 

We use IS A JET 7.80 [31] to calculate SUSY mass spectrum and decay modes. Models 

are required to have the cross section larger than 0.1 pb in the 14 TeV pp-collision so 

that the large enough number of events will occur at the LHC. The cross sections of the 

SUSY events are calculated by HERWIG 6.510 [35], [36]. We also require that the stau, 

winos and colored superparticles (i.e., sup, sdown and gluino) are heavier than 100 GeV, 

260 GeV and 450 GeV, respectively. In addition, we impose a vacuum stability condition, 

\A T \ tan /3 sin 28 T < 400 [37], where A T and 9 T are A-parameter and the mixing angle of 

f . For the scenario (A) and (B), f\ is the LSP in the MSSM-sector and other superparticle 

masses are larger than + 1 GeV. For the scenario (C), the same mass condition is 

imposed to superparticles except for Xv Xi 1S required to be the LSP while fi is the 

1 Since we pay particular attention to the case with small RPV couplings, RPV effects on the SUSY 
mass spectrum are negligible. 



4 



NLSP and their masses are degenerated within 1 GeV. Then, 2000 parameter points are 
obtained for each scenario. After the parameter generation, SUSY events are generated 
by HERWIG, and then they are passed through the fast detector simulation package 
PGS 4 [38] with a slight modification to treat stable stau. f\ is tagged as fi-like if its 
velocity is smaller than 0.9c, while f\ with high-velocity is regarded as //-like. 

As we discussed above, a difference of underlying scenarios will be imprinted in a 
missing energy distribution and a /i-like track distribution associated to the event with 
one fi-like track. Here, we consider two quantities which are useful to extract a nature 
of the underlying model. The first one is the difference between the missing p? averaged 
over the events with 1 fi-like track and that averaged over the events with 2 fi-like tracks, 

A(#r) = {fa\l fi-like track) - ($ T | 2 fx -like tracks). (4) 

The second one is the average of the number of large-pr muon-like tracks and fi-like 
tracks, 

(A^fi) = (the number of \x— like tracks with large— px and fi— like tracks). (5) 

We use events with at least one fi-like track. Then we count up the number of the 
fi-like tracks and the /i-like tracks with p? > m fl o □ Notice that these quantities well 
describe the underlying model when the staus produced from heavy superparticle cascades. 
Therefore we use the event in which the leading jet has px > 80 GeV. In addition, we 
require that at least one of the following conditions is satisfied: (1) a second jet with 
Pt > 60 GeV, (2) an electron with p T > 30 GeV, (3) a r-tagged jet with p T > 40 GeV. 

In Fig. [H we show the distribution of these quantities in three scenarios. In sce- 
nario (A) and (B), two staus arise in the final state. Even if we cannot distinguish a 
high-velocity stau from muons, the resulting /x-like track would have large px- Hence we 
expect (Nnfa) ~ 2. In these cases, sources of the missing energy are mainly neutrinos and 
miss-measurement of jet energies. If a correlation between the number of fi-like tracks 
and the missing energy sources is weak, we expect A(^r) ~ 0. On the other hand, the 



2 The mass of f% is directly measurable using its track information [2j [39] . Hence we expect that the 
mass determination is accurate enough by using a small number of the tracks [5] . 



The result docs not change drastically as long as we adopt the pt > 0(nif 
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Figure 1: Distribution of A(j/>t) and (N^fi) i n three scenarios. The stau fraction is at most 
Vo in the scenario (C). On the contrary, the fraction is 100% in the other scenarios. 



final state would contain not only ?x but Xi m scenario (C). Therefore, we expect that 
(-W^fi) < 2 and A (fa) > in this scenario. These expectations are indeed true as one can 
see from the figure, and the scenario (C) would be discerned from the other two scenarios. 

For the discrimination of the scenario (C) from the other scenarios, roughly speaking, 
we need to determine the (N^) and A (fa) with the accuracies of 0.2 and 50 GeV, 
respectively. We estimate the sizes of the statistical fluctuations. With 100 signal events, 
5(A($t)) — ±50 GeV and 5({N^ fl )) ~ ±0.1. In the figure, we indicate a fraction of the 
stau appeared in the final state at the parton level. Both A(^>t) and (N^) are strongly 
correlated to the stau fraction. It will be possible to determine the stau fraction when 
the large number of signal events are observed. 



4 Decay of Stopping stau at the Hadron Calorimeter 

At the LHC, the stau with small velocity (/?7 < 0.45) lose their kinetic energy by ionization 
and stops at the hadron calorimeter [10] . The ATLAS detector can identify the stopping- 
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stau events and can observe the energy deposit of hadrons produced by the decay of 
such stopping staus with a wide range of the lifetime 0(0.1) — O(10 10 ) sec at the hadron 
calorimeter [10]. Non hadronic objects from the stau are difficult to be observed. 

In the previous section, we found that by observing the LHC event topology, we can 
discriminate the degenerated mass scenario from the other two scenarios: (A) the gravitino 
LSP with fi-NLSP and (B) fi-LSP with RPV. In this section, we show that by studying 
the decay of stopping staus at the hadronic calorimeter, we can distinguish between the 
scenario (A) and the scenario (B) which are not distinguished by the LHC event topology. 

We use following two observable quantities in order to distinguish the two scenarios. 
The first one is P| iad defined by 



ihad the number of events where hadrons from the stau decay are observed 



(6) 



the number of the stopping staus 
The second one is the energy distribution of hadrons produced by the stau decay. Here we 



of the staus decay within the time that the detector is prepared to observe 
We also assume that we can observe all of the energy deposit of hadrons 



assume that al 
the stau decay 3 

produced from the stau decayl^l Under these assumptions, P^ ad is equal to the branching 
fraction of the stau into at least one hadron. We use the programs Madgraph4.0 jlQ] and 
TAUOLA2.6 [H] to simulate the energy distribution. 

In the case of (A), a stau decays into a gravitino as fi — > tG and the tau immedi- 
ately decays leptonically or hadronically. Because the event with no hadrons can not be 
triggered at the detector, we can observe only the events where the tau decays hadroni- 
cally. Therefore, P| iad in this scenario is equal to the hadronic branching fraction of tau 
P T (~ 0.65). The energy distribution of hadrons is shown by the black line in Fig. HJ 
In the figure, we assume that the stau is purely right handed, and the stau mass is 
set to be 250 GeV. We also show the energy distribution with the detector resolution 
AE/E = 150%/v/P/GeV in Fig. El 

On the other hand, in the case of (B), many types of the decay modes are possible 
depending on the coupling constant of R-parity violating term. Here we concentrate on 



4 For detailed discussion about how to observe the decay of the stopping stau with various lifetime, 
see [10]. 

5 We include all of the energy of quarks into the energy of hadrons though quarks may radiate some 
non-hadronic particles. The energy of ir° produced by tau-decay is also included in the energy of hadrons. 
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Figure 2: Energy distribution of hadrons in each scenario for = 250 GeV. The 
normalization is arbitrary. Mode 1 and Mode 3-6 are the cases with RPV decay. (See 
Tabled) 

the trilinear R-parity violating terms in the superpotential, 



W 7 ^ = -^KjkLiLjE k + X' ijk LiQjD k + -\" jk UiDjD k , (7) 

where Li and Qi are superfield of SU(2) doublet of left-handed leptons and left handed 
quarks respectively. E{, Ui and Di denote the right-handed superfield of leptons, up-type 
quarks and down-type quarks respectively. The decay pattern of the stau by the trilinear 
R-parity violating coupling has been already discussed in the other context [12]. We 
summarize the decay modes and the hadronic branching ratio of staus when one R-parity 
violating coupling dominates in Table [TJ We show the energy distributions corresponding 
to each modes in Fig. [2] where we assume that the stau is purely right handed. The 
smeared energy distributions are shown in Fig. El 

We find that we can distinguish the scenario (A) from Mode 2-6 by observing P^ ad . 
When O(10) of stau decays are observed, we can discriminate scenario (A) from Mode 
2-6 in 3-a confidence level. It is also found that we can discriminate the scenario (A) 
from Mode 1 by observing the energy distributions of hadrons. We check that by using 
the maximum likelihood analysis to the smeared energy distributions of scenario (A) 
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Figure 3: Same as Fig J2] but with the energy resolution AE/E = 150%/ \/E/GeV. 



Mode 


coupling 


decay modes 


hadronic branch 


Mode 1 


^121 = — A211 
^122 = — ^212 


r + efeji7 k 


Pr 




^131 = — A311 








Al32 = — A312 


4"j 




Mode 2 


^231 = — A321 
A232 = — A322 
^123 = — ^213 


4 v i 


< 1 


Mode 3 


^133 — — -^313 


e+v T 


P T /2 (purely right handed) 




-^233 — — -^323 


efi7 T 

T + h>i 

r + i>i 


I 

P T (purely left handed) 


Mode 4 


\l 

A llm 








A 2lm 


T + didji7 k 
T + didjU k 


1 


Mode 5 


\l 

A 3lm 


Ujd k 


1 


Mode 6 


A lmn 


T + Uidjd k 
T + Uidjd k 


1 


Scenario (A) 


T+G 


Pr 



Table 1: Decay mode and hadronic branching ratio of when one R-parity violating 
coupling dominates. We also show those of the scenario (A) for the reference. (ei,e2) = 
(e,/x), (vx,v 2 ) = (v e ,Vp,), {v,x,U2,u$) = (u,c,t) and (d 1 ,d 2 ,d 3 ) = (d,s,b). 
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and Mode 1 in scenario (B), the discrimination of the two energy distributions in 3-a 
confidence level is possible by observing O(10) of stau decays. Thus, by observing 0(10) 
of stau decays, we can distinguish the scenario (A) with (B) in 3-cr confidence level. It is 
also checked that even when the several modes contribute, we can distinguish the scenario 
(A) with (B) by using the hadronic branching fraction and the energy distributions. In 
the realistic detector system, only the hadrons with the energy more than 20 GeV can be 
observed, however the conclusion does not change by the energy cut. 

In this paper, we do not consider the bilinear term tiLiH u in RPV case where H u is the 
superfield of up-type Higgs. In the scenario with the bilinear couplings, there are many 
types of the decay mode such as fi — > tu, f\ — > tb and fi — > WZ. We can discriminate 
the two scenarios by using the same way in the large parameter region. In some cases, 
unfortunately, a stau mainly decays by the process f\ — > tu or fi — > ru, and therefore 
the hadronic branching ratio and the energy distribution of stau is almost the same as 
the stau NLSP case. In such cases, it is difficult to distinguish the scenario (A) from the 
scenario (B). 

The higher dimensional RPV operators are also possible. With the higher dimensional 
RPV operator, we expect that by comparing the energy distribution, it is possible to dis- 
criminate the two scenarios because it is expected that the decay using higher dimensional 
operator tend to produce more than 2 particles. The study of the case with the bilinear 
terms and the higher dimensional operators are for future work. 

5 Conclusion and Discussion 

In this paper, we have discussed how to prove the origin of the stau-longevity at the LHC 
experiments. There exist three distinct mechanisms which make the lighter stau long- 
lived: weak interaction between fi-NLSP and LSP, fx-LSP with small R-parity violation 
and suppressed fi decay into the LSP due to the mass degeneracy. 

If the masses of the fi-NLSP and x?-LSP are sufficiently degenerated, the stau does 
not promptly decay into the neutralino due to the kinematical suppression. In such a 
scenario, heavy superparticles produced at the LHC will decay into not only the stau but 
the neutralino. This feature can be probed via SUSY event topology. In particular, the 
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missing energy distribution and the number of high-p<r muons provide useful information. 
They have been investigated in section [31 and it is shown to be possible to check the 
degeneracy between the stau and the neutralino at the early stage of the LHC experiment. 

The decay of the stopping stau at the calorimeter provides another important infor- 
mation on the origin of the longevity as discussed in section HJ If the stau is the NLSP 
and there is the weakly interacting LSP, e.g., gravitino or axino, f\ will produce r in its 
decay into the neutral LSP. In the case of the stau-LSP with .R-parity violation, there 
are various stau decay modes depending on the RPV couplings but they generally posses 
different final states from the fi-NLSP case. A hadronic branching fraction of the stau 
decay can be used to discriminate the fi-NLSP case from some of RPV cases. Also, an 
energy distribution of hadrons from stau decay is useful. Using these informations, we 
could discriminate the cases where RPV is induced by the trilinear coupling from the 
fi-NLSP scenarios. 

In this paper, we assume that all the decays of trapped staus occur in operation time of 
the detector, and that the detection efficiency for the hadronic decay of the stau trapped 
in the calorimeter is 100 %. In practice, there need be some corrections to estimate the 
hadronic branching fraction because some of the decays occur in the dead time of the 
detector, depending on the lifetime of the stau. Although it would be also necessary 
to consider the detection inefficiency, it is expected that such effects do not change the 
present discussion drastically. Detailed treatments of such effects are beyond the scope of 
this paper. 

Let us comment on the application range for the present method. The method of the 
event topology requires 0(100) events in order to discriminate the models. The method 
of the measurement of the stopped stau decay in the calorimeter needs O(10) stopping 
events. Such stopping events are typically 0(1) % of all SUSY events [10J. Therefore, 
the present method is applicable to the case that < 500 GeV or rn^q < 1500 GeV 
for the integrated luminosity 10 fb _1 at y/s = 14 TeV at the LHC. 
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